Abstract. The metal line profiles of different ions ob- 
served in high Hi column density systems [-/V(Hi) > 10^^ 
cm^^] in quasar spectra can be used to constrain the ion- 
ization structure and kinematic characteristics of the ab- 
sorbers. For these purposes, a modified Monte Carlo In- 
version (MCI) procedure was applied to the study of three 
absorption systems in the spectrum of the HDF-South 
quasar J2233-606 obtained with the UVES spectrograph 
at the VLT/Kueyen telescope. The MCI does not confirm 
variations of metal abundances within separate systems 
^vq " which were discussed in the literature. Instead, we found 
that an assumption of a homogeneous metal content and a 
unique photoionizing background is sufficient to describe 
■ the observed complex metal profiles. It was also found that 
P5 I the linear size L and the line-of-sight velocity dispersion 
, ^ ■ tTv measured within the absorbers obey a scaling relation, 
namely, CTv increases with increasing L, and that metal 
abundance is inversely proportional to the linear size of 
the system: the highest metallicity was measured in the 
system with the smallest L. 



00 



o 

O 

:^ 

o 



X 
J3 



Key words: Cosmology: observations - Line: formation 
- Line: profiles - Galaxies: abundances - Quasars: absorp- 
tion lines - Quasars: individual: J2233-606 



A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 

missing; you have not inserted them 



ASTRONOMY 

AND 
ASTROPHYSICS 



Metal abundances and kinematics of quasar absorbers 

I. Absorption systems toward J2233— 606* 

S. A. Levshakov^**, I. I. Agafonova^, M. Centurion^, and I. E. Mazets^ 

^ Division of Tlieoreticai Astrophysics, National Astronomical Observatory, Mitaka, Tokyo 181-8588, Japan 
^ Department of Theoretical Astrophysics, loffe Physico- Technical Institute, 194021 St. Petersburg, Russia 
^ Osservatorio Astronomico di Trieste, Via G. B. Tiepolo 11, 34131 Trieste, Italy 



Received 00 / Accepted 00 



1. Introduction 

Absorption systems in quasar spectra provide unique in- 
formation on the intervening intergalactic matter (IGM) 
up to redshift 2: ~ 6, back to the time when the Universe 
was less than 7% of its present age. High resolution spec- 
troscopic observations available nowadays at large tele- 
scopes open new opportunities to investigate the physical 
nature of quasar absorbers. Reliable data on the chemical 
composition of the IGM and on the physical characteris- 
tics (like velocity and density distributions, volumetric gas 
density, kinetic temperature, ionization structure etc.) of 
the absorbers is an important clue to our understanding 
of galaxy formation, chemical evolution of the IGM and 
the origin of the large-scale structure. 

In recent investigations much attention find the metal 
systems which are the absorbers exhibiting as a rule nu- 
merous lines of low (like Hi, Cii, Sin, Mgii, Alil) and 
high (like Cm, Niii, Sim, Civ, Siiv, Nv) ionized species. 

Presence of metals provides a unique opportunity to 
study the physical conditions of matter at early epochs. 
Unfortunately, the computational methods usually ap- 
plied to high resolution spectra lie quite often behind the 
quality of observational data and fail to extract from them 
all encoded information. The common processing method 
consists of the deconvolution of complex absorption pro- 
files into an arbitrary number of separate components (as- 
suming a constant gas density within each of them) which 
are then fitted to Voigt profiles. However, in many cases 
this procedure may not correspond to real physical condi- 
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tions: observed complexity and non-similarity of the pro- 
file shapes of different ions indicate that these systems are 
in general absorbers with highly fluctuating density and 
velocity fields tightly correlated with each other. Too high 
or too low gas temperatures, extremely varying metallici- 
ties between subcomponents, exotic UV background spec- 
tra and other physically badly founded outcomings may 
be artifacts arising from the inconsistency of the applied 
methods (see examples in Levshakov et al. 1999; Lev- 
shakov et al. 2000b, hereafter Paper I). 

In Paper I we developed a new method for the QSO 
spectra inversion, - the Monte Carlo Inversion (MCI), - 
assuming that the absorbing region is a cloud with uniform 
metallicity but with fluctuating density and velocity flelds 
inside it. This computational procedure which is based on 
stochastic optimization allows us to recover both the un- 
derlying hydrodynamical fields and the physical parame- 
ters of the gas. First application of the MCI to the analysis 
of the Zabs = 3.514 system toward Q08279+5255 (Lev- 
shakov et al. 2000a) has shown that the proposed method 
is very promising especially in the inversion of complex 
absorption spectra with many metal lines. 

In this paper we start a new comprehensive survey 
of the metal systems for which high resolution and high 
signal-to-noise spectra are available. We present here the 
results for three absorption systems {zahs = 1.87, 1.92 and 
1.94) from the spectrum of the quasar J2233-606 which 
have been already studied by Prochaska & Buries (1999), 
and D'Odorico & Petitjean (2001, hereafter DP) using the 
common Voigt fitting method. We re-calculate these sys- 
tems using the MCI in order to compare the applicability 
of both approaches and to show up their restrictions. 

The structure of the paper is as follows. In Sect. 2 the 
data sets used in the MCI analysis are described. Sect. 3 
contains the details of the applied computational proce- 
dure. The results obtained for each of the mentioned above 
systems are presented in Sect. 4. Conclusions are reported 
in Sect. 5. In Appendix the general equation of the en- 
tropy production rate is given which is used to calculate 
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Fig. 1. Hydrogen and metal absorption lines associated with the Zahs — 1-92 system toward J2233-606 (normalized intensities 
are shown by dots with la error bars). The zero radial velocity is fixed at 2 = 1.92595. Smooth lines arc the synthetic spectra 
convolved with the corresponding point-spread functions (FWHMvlt = 6.7 km s"'^, FWHMhst = 10.0 km s~^) and computed 
with the physical parameters from Table 1. Bold horizontal lines mark pixels included in the optimization procedure. The 
normalized Xmin = 1-10 (the number of degrees of freedom u = 1450) 

the density and velocity configurations along the line of AA 2275 — 3118 A) obtained with the HST/STIS (Savaglio 
sight exhibiting minimum dissipation. 1998). 



2. Data 

High-quality data of the Hubble Deep Field South (HDF- 
S) quasar J2233-606 (zem = 2.238, B ~ 17.5) were ob- 
tained during the commissioning of the UVES on the 
VLT 8.2m Kucycn telescope at Paranal (Chile) in Oc- 
tober 1999. The resolving powers at which the spectra 
were recorded in the spectral ranges AA = 3050 — 10000 A 
( J2233-606) and AA = 3700 - 5000, 6670 - 9100 A (QOOOO- 
2620) were R ~ 45000 and R 49000, which correspond 
to velocity resolutions of FWHM ~ 6.7 km s~^ and ~ 6.1 
km s"'^, respectively. 

The data reduction and the identification of metal 
absorption-line systems in the J2233-606 spectrum are 
reported in Cristiani & D'Odorico (2000). In our study 
we also used the J2233-606 echelle spectrum {R ~ 30000, 



3. Computational procedure 

The complete description of our computational procedure 
is given in Paper I. Here we summarize briefly basic model 
assumptions and emphasize new details recently included 
in the MCI. 

We assume that all lines observed in a metal system 
arise in a continuous absorbing gas slab of a thickness 
L (presumably the outer region of a foreground distant 
galaxy). The absorber exhibits a fluctuating gas density 
and a mixture of bulk motions such as infall and out- 
flows, tidal flows etc., resulting in a stochastic velocity 
field. Metal abundances are assumed to be constant within 
the absorber and gas is supposed to be optically thin for 
the ionizing UV radiation. 



4 



S. A. Levshakov et al.: Metal abundances and kinematics of quasar absorbers 



Within the absorbing region the radial velocity v{s) 
and total hydrogen density nu{s) distributions along the 
line of sight are the same for all ions. In the computa- 
tional procedure these two random fields are represented 
by their sampled values at equally spaced intervals Ax {x 
is dimensionless radial coordinate s/L), i.e. by the vec- 
tors {vi, . . . ,Vk} and {ni, . . . , n^} with k large enough to 
describe the narrowest components of complex spectral 
lines. 

Further wc suppose the thermodynamic and ioniza- 
tion equilibrium at each computational point along the 
sightline which means that fractional ionizations of dif- 
ferent ions are determined exclusively by the gas density 
and vary from point to point. These fractional ioniza- 
tion variations are just the cause of the observed diver- 
sity of profile shapes. To calculate the kinetic temperature 
and fractional ionization of ions the photoionization code 
CLOUDY (Ferland 1997) was used. 

The inputs to CLOUDY are the dimensionless ion- 
ization parameter U = Uph/nn (riph - the number den- 
sity of photons with energies above 1 Ry), metallicity and 
the background ionizing spectrum for which the Haardt- 
Madau spectrum (HM) was adopted (Haardt & Madau 
1996). The number density of the ionizing photons for this 
spectrum 

nph = -r'^912 / ^— — , (1) 
ch J^^ \J912J y 

equals at z = 2 to 2.26 x 10~^ cm~^ [here c, /i, and 

are the speed of light, the Planck constant, the frequency 

of the hydrogen Lyman edge, and the specific intensity (in 

ergs cm~^ s^^ sr^^ Hz^-'^) at frequency v, respectively]. 

Fractional ionization curves T([/) were computed 
with CLOUDY for solar abundance pattern and different 
metallicities and then included in the MCI code to calcu- 
late the optical depths for ions involved in the fitting. If 
the obtained solution revealed the abundance pattern dif- 
ferent from the solar one, T(?7) were recalculated for this 
new pattern and all computations repeated. It should be 
noted, however, that differences of 0.2 -=-0.3 dex from solar 
values infiuence the fractional ionizations only weakly. 

The values of velocity and density at subsequent com- 
putational points are considered to be correlated and are 
described by means of Markovian processes. In particular, 
the velocity is computed as follows: 
v{x Aa;) = /v v{x) -\- e{x Aa;) , (2) 
where /v = i?v(Ax)/cry, Ry being the correlation between 
the velocity values at points separated by a distance Ax, 
i.e. J^v = kp^"^ + Ax^v{x)), Gy the velocity dispersion of 
bulk motions, and e a random normal variable with zero 
mean and dispersion 

o-e.v = o-v^l - /v • (3) 
Density is supposed to have a log-normal distribution 
(characterized by the mean tiq and the second central mo- 
ment (Ty of dimensionless variable y = nu/rio ) and is de- 
scribed with help of an auxiliary Markovian process ^{x): 
v{x + Ax) = v(x) + e(x + Ax) , (4) 
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Fig. 2. Computed velocity (panel a) and density of gas (panel 
b) and ions (panels c-j) distributions along the line of sight 
for the system at 2 = 1.92595 toward J2233-606. Shown are 
patterns rearranged according to the principle of minimum en- 
tropy production rate (see text) 



where fi^ = Ri,{Ax)/al, Ri, being the correlation between 
the values of v at points separated by a distance Aa;, 

the logarithmic density dispersion, Ov 



'ln(l + aj) and 
e a random normal variable with zero mean and dispersion 

Having defined v{x), the total hydrogen density can 
be obtained as 

nH(ar) =noe^(^)-i'^' . (5) 
The optical depth Ta,i(A) at wavelength A for element 
'a' in ith ionizing stage is calculated from the equation: 

Ta,,(A) = fco^aiVoe-i-' 



e^(^) Ta,i[C/(a;)] ^x\A\t,{x),v{x)\ da; , 



(6) 



where fco is a constant, iVo = i^q L - the expectation vahic 
of the total hydrogen column density, Za = ng^/nu - the 
relative abundance of element 'a', Ta,i[C^(^)] - the frac- 
tional ionization for stage i, $a ~ the profile function, 
AAd - the Doppler width of the line. 
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Fig. 3. Density-weighted velocity distribution functions, p(v), 
for Hi, Cii, Sin, Niii, Sim, Civ, and Siiv as restored by the 
MCI procedure in the z = 1.92595 system toward J2233-606 



Before being compared with the observed spectrum, 

the synthetic intensities cxp[— t(A)] are convolved with the 
spectrograph point-spread function. 

Thus the proposed model is fully defined by specify- 
ing the following values: the velocity vector {vi}, the total 
hydrogen density vector {rij}, the total hydrogen column 
density A^O; the mean ionization parameter Uq, the radial 
velocity dispersion (jy , the density second central moment 
fTy, the element abundances Za and the correlation coeffi- 
cients /v and /i/. The common least-squares minimization 
(LSM) of the objective function is used to estimate the 
model parameters'^. 

The computational procedure itself consists of 
two stcips: firstly a point in the parameter space 
(A^O) Uo, CTv, CTy , Za) is chosen and then an optimal config- 
uration of {vi} and {ui} for this parameter set is searched 



* The LSM is computationally identical to the standard x^- 
minimization, but it can be used to evaluate the uncertainties 
of the fitting parameters even in cases of correlated measure- 
ments. One must realize that the rebinning used in the data re- 
duction procedure introduces a correlation between data point 
values (Levshakov et al. 2002) 



for. Correlation coefficients arc considered as external pa- 
rameters and remain fixed during the calculations. 

The optimization of {vi} and {rii} is the most time- 
consuming part of the procedure and needs an effective 
algorithm to achieve a quick and stable convergence of the 
computations. In the MCI we use the simulated anneal- 
ing with Tsallis acceptance rule (Xiang et al., 1997) and 
an adaptive annealing temperature choice. Namely, the 
annealing temperature Ta at iteration -|- 1 is decreased 
according to following equation: 

rpk+1 £a /yv 

^l + P^r' ^' 

where AEk/Ek is the relative energy (= the sum of 
squares of the deviations, S) variation on step fc, r is the 
acceptance rate, i.e. the ratio of the accepted trials to the 
total trial number, and /? is a constant of order 1. 

The calculation of the uncertainty ranges for the fit- 
ting parameters is in our method not so straightforward 
as a simple inversion of the Hesse matrix since the veloc- 
ity and density distributions represent additional degrees 
of freedom and widen, in general, the confidence intervals. 
However, these {vi} and {rii} distributions themselves are 
nuisance parameters and should be 'integrated out' when 
one computes the errors for the other parameters. To es- 
timate the confidence levels, the following procedure can 
be applied: the values of the physical parameters in the 
vicinity of the global minimum of the objective function 
are chosen at random and then the optimal density and 
velocity distributions are computed. Assuming that the 
probability of each parameter set can be linked to the de- 
rived S value [e.g. as exp(— 5)] we can estimate from the 
obtained sample the joint probability density function for 
parameters and hence calculate all necessary statistical 
moments. 

It should be noted, however, that the reliable estima- 
tion of this multidimensional function requires a very large 
sample which is quite time consuming. In our case the 
exact estimation of confidence levels is not very crucial 
taking into account the intrinsic uncertainties in atomic 
data (e.g. Savin 2000) or unknown shape and intensity of 
the local background ionizing radiation. Because of this 
we restricted our samples to a few dozens of points and 
estimated the accuracy of the fitting parameters only ap- 
proximately. 

The recovered density and velocity patterns arc not 
unique - many configurations are possible with compa- 
rable probability. But all these configurations have the 
same density-weighted velocity distributions which actu- 
ally determine the observed line shapes (see Paper I). As 
already mentioned above, we represent these random fields 
by their values sampled at equally spaced intervals Ax. 
In order to compare the calculated patterns we rearrange 
these values in such a way that the final configuration 
exhibits a lowest rate of entropy production: according to 
the Prigogine theorem (Prigogine 1967), this configuration 
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Fig. 4. Same as Fig. 1 but for the Zs,ha ~ 1.94 system toward J2233-606. The zero radial velocity is fixed at 2 = 1.942616. 
The corresponding physical parameters are listed in Table 1. The normalized Xmin = 1-08 (the number of degrees of freedom 
z/ = 888) 



has the minimal dissipation and, hence, is more stable and 
more probable. All necessary equations to calculate the 
entropy production are presented in Appendix. We stress, 
however, that configurations produced on the base of such 
rearrangement should in no case be considered as some- 
thing final - they represent only the (most) probable case 
of the density and velocity distributions along the line of 
sight and are used here exclusively for illustrations. 

4. Results on individual systems 

4.1. J2233-606, ^;abs ^ 1-92 

The system at ^;abs — 1-92 toward the quasar J2233-606 
has saturated hydrogen lines (from Ly-a up to Ly-8) and 
metal lines of Cii, Sin, Sim, Civ, and Siiv. The results 
obtained with the MCI are presented in Table 1 and illus- 
trated in Figs. 1 and 2. Parts of profiles included in the 
minimization are marked by horizontal lines at the panel 
bottoms in Fig. 1. Profiles of the doublet Nv A1238 A and 
Nv A1242 A were calculated later using the obtained ve- 
locity and density distribution and the metallicity derived 
from the fitting of N III A989 A. It is seen from Fig. 1 that 
most spectral features can be well represented assuming 
uniform metallicities and a common HM UV background. 



Fig. 2 demonstrates the distribution of the radial velocity 
and gas density (panels a and b) along the line of sight 
(rearranged in accord with the principle of minimal en- 
tropy production rate). The density distributions for the 
ions involved in the optimization are shown in panels c-j, 
whereas the density-weighted velocity distributions which 
determine the shapes of the spectral lines are presented 
in Fig. 3. This figure shows that the density-weighted ve- 
locity distributions for low ions Cii and Sin are similar, 
but differ from those for high ions Civ and Siiv. These 
distributions easily explain why the lines of C li and Si li 
look very much alike and why their centers are displaced 
by At; ~ 4 km s~^ (DP) with respect to Civ and Si IV. 

The study of this system by DP. who used the stan- 
dard Voigt profile fitting, produced comparable column 
densities (albeit 20-50% smaller). However, the metallici- 
ties obtained by DP for two main clouds (at w = km 
and V = —43 km s~^) differ nearly by two orders of mag- 
nitude: [X/H] = -0.9 and -2.7, respectively (in our case 
[X/H] ~ —2.2 for the whole system). As shown in Pa- 
per I, the Voigt fitting may in general yield correct col- 
umn densities when applied to unsaturated lines, but the 
mean U and, hence, the ionization corrections may not 
be unambiguous. Therefore, the conclusion made by DP 
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Fig. 5. Computed velocity (upper panel) and gas density 
(lower panel) distributions along the line of sight for the sys- 
tem at z = 1.942616 toward J2233-606. Shown are patterns 
rearranged according to the principle of minimum entropy pro- 
duction rate (see text) 

that the z^hs — 1-92 system contains 'a region of intense 
star-formation activity' may not be well justified since this 
result is model dependent. 

The values of the average gas density no and kinetic 
temperature Tkin, and the cloud thickness L estimated 
in our model (sec Table 1) are typical for the Ly-a sys- 
tems discussed in the literature (e.g., Giallongo & Petit- 
jean 1994; Viegas et al. 1999; Prochaska & Buries 1999; 
Chen ct al., 1998; Chen et al. 2001). Low mctallicity for 
the whole system ([X/H] < —2.0) and its dimension of 
20 kpc imply that this system can originate in a galactic 
halo or in a large scale structure object. 

12. J2233-606, Zabs ^ 1-94 

This system exhibits a plenty of metal lines in different 
ionization stages. The metal profiles are not very complex 
and extend over the velocity range from —100 km s^"'^ to 
100 km s~^. Results obtained with the MCI are presented 
in Table 1 and shown in Figs. 4 and 5. As in the previous 
system, most absorption features can be well described 
with uniform metallicities and a common HM spectrum. 
The Ly-a profile is contaminated by the forest absorp- 
tion in the blue and red wings and therefore the Ly-a 
absorption feature was not involved in the analysis. The 
profiles of MgiiA2796 A and AliiA1670 A were com- 
puted later using the derived velocity and density distri- 
butions. MgiiA2796 A is contaminated by a telluric line 
and this explains the difference between the computed and 
observed profiles. The synthetic and observed profiles of 
AliiA1670 A show much more pronounced discrepancy. 
Fractional ionisation curves for Al II and Al iii were com- 
puted with CLOUDY. These curves allowed us to fit the 
Aim doublet quite well with the Al abundance similar 
to that obtained for the other metals. However, when the 
Alii profile was included in the fitting, the Al mctallic- 
ity differed by order of magnitude from the other metals. 



Besides it was impossible to fit adequately the Al ill dou- 
blet. Similar behaviour of Al was reported also by DP who 
noted that 'the recombination coefficients used to com- 
pute the aluminium ionisation equilibrium [in CLOUDY] 
are probably questionable'. 

Column densities derived by DP coincide well (within 
15%) with that obtained in our procedure except for the 
saturated SiiiiA1206 A line for which the Voigt fitting 
gave nearly 2 times lower value. The abundances esti- 
mated in DP scatter again from component to component, 
but nevertheless they conclude that 'the gas in this system 
is likely of quite high metallicity (larger than 0.1 solar)'. 

Similar to the Voigt fitting, the MCI also delivered for 
this system high metal abundancies; one third solar for 
carbon and silicon and nearly two times lower for nitro- 
gen, magnesium and aluminium. Taking into account this 
result and a compact dimension (~ 5 kpc, see Table 1) of 
the absorbing region we come to the same conclusion as 
Prochaska & Buries (1999) did: the system at z = 1.94 can 
hardly be a large scale structure object (like a filament or 
a wall) and should be related to a galactic system (may 
be a region of intense star formation). 

4.3. J2233-606, Zabs ^ 1-87 

This is the most interesting system from the family of the 
absorbers at 2; = 1.9 toward J2233-606. The metal line 
profiles show a rather complex structure extending over 
the velocity range of about 700 km s~^. Some of these pro- 
files are severely blended that hampers the unique Voigt 
profile deconvolution (e.g. DP assumed 17 components to 
describe metal profiles). 

The MCI code turned out to be much more robust and 
was able to recover the self-consistent line profiles even 
under such unfavourable conditions. The physical param- 
eters which the MCI delivered for the z = 1.87 system 
together with the underlying velocity and density distri- 
butions are presented in Table 1 and in Figs. 6 and 7. It 
is seen from Fig. 6 that like in the previous two systems 
all lines are well described with a single parameter set, 
uniform metallicities and a common HM UV background. 
The blue wing of the Ly-a line is contaminated by the 
forest absorption as is clearly seen from the Ly-/3 and Ly- 
7 profiles. The synthetic profile of the O VI A 1031 A line 
was calculated later using the derived best fitting param- 
eters and the oxygen abundance [0/H] = — 1.0 (which is 
about 3 times over the other element abundances from 
this system). Even with the increased abundance the syn- 
thetic profile of O vi is still much weaker than the observed 
intensities. This discrepancy rules out the ionization of 
Ovi by the adopted background radiation. Taking into 
account that all other elements have been well described 
with a given HM spectrum and that the coUisional ion- 
ization of oxygen can hardly be effective at low densities 
(riH ~ 10~^) and temperatures of ^ 25000 K, this result 
seems to favor the interpretation that the Ovi ion and 
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Fig. 6. Same as Fig. 1 but for the ^abs — 1-87 system toward J2233-606. The zero radial velocity is fixed aX z = 1.87008. 
The corresponding physical parameters are listed in Table 1. The normalized Xmin = 1-60 (the number of degrees of freedom 
z/ = 1158) 



the other ions do not arise in the same gas (Kirkman & 
Tytlcr 1999; Reimers et al. 2001). 

According to our results, the absorber at ^ = 1.87 
could be a large size cloud with very high velocity dis- 
persion. Its estimated linear size of 80 kpc is consistent 
with dimensions of extended gaseous envelopes observed 
around galaxies at 2; < 1. In these envelopes, Mgii ab- 
sorption is the dominant observational signature at the 
distancies up to a few tens of kiloparsecs (Bergeron & 
Boisse 1991), whereas highly ionizied species like CiV are 
observed at distances of at least 100 kpc from galactic 
centers (Chen et al. 2001). Since the extended structure 
of the same order of magnitude is observed at 2 = 1.87, 
we may conclude that this system arises in the external 
halo at large galacto centric distances. 

5. Discussion and conclusions 

A main goal of this work was to investigate the reliabil- 
ity of the physical parameters and dynamical characteris- 
tics of the metal absorption systems obtained by means of 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Space Coordinate, x — s/L 

Fig. 7. Same as Fig. 5 but for the ^abs — 1-87 system toward 
J2233-606 

the standard Voigt fitting procedure and by the modified 
Monte Carlo Inversion algorithm. For comparison we used 
the recently obtained results on the Voigt fitting analysis 
of three systems toward J2233-606 (D'Odorico & Petit- 
jean 2001). 
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Wc found that both approaches dchvcr similar total 
column densities of unsaturated metal lines. The saturated 
profiles may, however, be treated differently (e.g., the Sim 
A 1206 A line in the 2; = 1.94 system). 

We also found that metal abundances based on the 
Voigt deconvolution procedure differ considerably from 
those obtained by the MCI. For instance, instead of fluc- 
tuating metallicitics found in the absorbers toward ,12233 
606 by DP, the MCI shows that an assumption of a ho- 
mogeneous metal abundance for the whole system under 
study is quite sufficient to represent all observed features. 

New and principal results which can be obtained only 
with the MCI procedure arc the kinematic characteristics 
of the absorbers. We estimated selfconsistently for the first 
time the density and velocity dispersions along the sight- 
lincs within the absorbers and calculated the total hy- 
drogen column and volumetric densities which gave us a 
direct measure of their linear sizes. The found dimensions 
of _L ~ 5 kpc to i ~ 80 kpc arc in good agreement with 
measurements of the extended gaseous envelopes around 
the nearby galaxes which were probed by the Mgll ab- 
sorption lines (Bergeron & Boisse 1991) and by the CiV 
lines (Chen et al. 2001). 

A new issue obtained in the MCI analysis is a scaling 
relation. Namely, we found that the linear size L shows a 
positive correlation with the line-of-sight velocity disper- 
sion (Tv, i.e. the higher L, the larger is observed (sec 
Table 1). Although our sample is still too small to carry 
out statistical analysis of this correlation, the scaling ten- 
dency is of the same kind that can be expected for viri- 
alized systems. The velocity dispersion is closely related 
to the total mass of the system in a stationary state (cf. 
the scaling law known as the fundamental plane for ellip- 
tical galaxies). Taking into account that the scaling laws 
are different for different types of objects (see, e.g., Fig. 2 
in Mallen-Ornclas et al. 1999), future statistical analysis 
may allow us to classify absorbers at different redshifts. 

It is also interesting to note another scaling relation: 
we observe systematically higher metal abundance with 
decreasing L, and vice verse, the higher L, the lower metal- 
licity is deduced. If L reflects the linear size of a distant 
absorber, then we may conclude that a compact absorber 
has, presumably, higher metal content as compared with 
an extended one. 
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Appendix A: Minimum entropy production rate 



General equation for the entropy production rate is given 
by (Landau & Lifshits 1987): 
dS_ _ f K(gradT)2 
'dt ~ 



J^2 

V / dvi^ dvk_ 
2T \dxk dxi 

|(divv)2dl/. 



AV- 



3 ox 



(A.1) 



Here dS/dt is the rate of the total entropy production in 
the system, the velocity v and temperature T are functions 
of the coordinate (.x,y, z). the volume element is denoted 
by dV , the thermal conductivity, dynamic viscosity, and 
the second viscosity coefficients - by k, r?, and ^, respec- 
tively. 

The terms included in this equation account for only 
hydro- and thermodynamic processes since we assume the 
ionization balance in each point of the region which means 
that the radiative heating and cooling do not contribute to 
the entropy production. The second viscosity equals zero 
for dilute monoatomic gases (Chapman & Cowling 1970) 
so we omit it from further consideration. 

We also assume that heat conductivity k and viscosity 
rj are dominated by turbulence, and hence the Prandtl 
number is about unity (Monin & Yaglom 1975): 



Pr 



1 



(A.2) 
1% 

where Cp is a specific heat capacity at constant pressure. 
The relation (A.2) is valid for both the ionized and neutral 
gas. 

From observations, only one component Vx (along the 
sightline) of the velocity vector v is known. Therefore 
we are compelled to neglect in (A.l) all terms includ- 
ing derivatives other than dvx/dx. The second right hand 
term in (A.l) can be re- written in the form: 



r] 


f dvi 


l_ dvk 


2T 


\dxk 


dxi 


11 


( dv. 


l_ dvk 


2T 


\dxk 


dxi 



6 oxj 



dV 



dV 

IL 

2T 



(divv)^ dV (A.3) 



After all these assumptions, we obtain the following 
simplified ID form for the entropy production rate 
dS_ _ 
'di ~ 

2 




+ Pr 



A 



CpTo T 



dvx 
dx 



dx 



(A.4) 



where Vx = Vx/(Jv, T = T/Tq, with and Tq being scales 
(characteristic values) for velocity and temperature inside 
the considered region, ^ is a constant of about 1, and S 
is the entropy per unit area. 

For a monatomic ideal gas Cp = 5nR/2 in is the num- 
ber of gram- moles and R is the universal gas constant). 
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Table 1. Physical parameters of the metal absorption systems toward J2233-606 derived by the MCI procedure 



raramet er 


Zahs — i.y4 


y . ~ 1 Q9 

^abs — 


Zgbs — '-■01 


Mean ionization parameter, Uo 


8.79-2 (~ 15%)* 


2.99-2 (~ 20%) 


0.12 (~ 10%) 


Total H column density, A'h, cm~^ 


1.41i9 (~ 15%) 


2.6920 (~ 10%) 


1.2320 15%) 


Velocity dispersion, av, km s~^ 


35 (~ 15%) 


104 (~ 20%) 


179 (~ 10%) 


Density dispersion, o"y 


1.50 (~ 15%) 


2.31 (~ 10%) 


1.31 (~ 10%) 


Chemical abundances^: 








Zc 


9.8-b(< 10%) 


2.0-6 {■^ 20%) 


8.7-6 (< 10%) 


Zn 


1.1-5 {-^ 15%) 


5.5-7 {■^ 20%) 




Zms 


6.8-6 (=i 10%) 




2.0-6 ( < 10%) 


Zai 


5.1-7 (=i 15%) 






Zsi 


1.0-5 ( < 10%) 


2.5-7 (=i 20%) 


1.7-6 (< 10%) 


[Zc] 


-0.53 


-2.21 


-1.58 


[Z^] 


-0.88 


-2.17 




[Zms] 


-0.75 




-1.27 


[Zm] 


-0.79 








-0.54 


-2.14 


-1.31 


Column densities, cm ^ : 








iV(Hi) 


(2.50 d= 0.05)i6 


(2.2±0.25)i7 


(2.65±0.28)i6 


iV(Cii) 


(7.49±0.11)i3 


(1.80 ±0.20)13 




Af(Mgii) 


(6.04±0.22)i2 




(1.96 ±0.21)12 


iV(Siii) 


(8.79 ± 0.25)i2 


(2.40±0.16)i2 


(2.56 ±0.20)12 


N(Ciu) 






(5.36 ± 0.25)i4 


Af(Nni) 


(1.62±0.09)i4 






iV(Aliii) 


(1.63±0.08)i2 






iV(Siiii) 


(9.33±0.51)i3 


(2.90±0.18)i3 


(4.94 ±0.25)13 


A''(Civ) 


(3.93 ± 0.16)i4 


(4.78 ± 0.31)13 


(3.66 ± 0.13)14 


A^'fSiivl 

J.T \ kJllV I 


^5 58 ± 35)15 


(8.50 ±0.50)12 


(4.13 ±0.11)13 


iV(Nv) 


< 3.54i2 






Hydrogen number density, no, cm~^ 


-1-3 


-5-3 


-5-4 


Mean kinetic temperature, K 


143 


253 


253 


Minimum kinetic temperature, K 


113 


I63 


I83 


Maximum kinetic temperature, K 


173 


4O3 


31.53 


Linear size, L, kpc 


~ 5 


~ 20 


~ 77 


Mass, M*, in Mq 


~ 1.46 


~ 4.58 


-2.69 


Mass, M*, in Mq 


~ 7.08 


~ 2.O10 


~ 3.O11 


* Values like 'oa,' mean a x 10"^; shown 


in parenthesis are 


uncertainties at the la c.l. 



tZx = X/H, [Zx] = log(Zx) - log(Zx)0 

■^Assuming spherical geometry; Using the virial theorem 



and if the gas is fully ionized then k oc T^/^ (e.g. Lang 
1999). 

Given the values of k, Pr, Tq, Cp, and Uv, we 
can rearrange the computational points {vi, . . . ,Vk} and 
in such a manner, that minimum of (A. 4) 
will be achieved. Minimization of (A. 4) was carried out 
by means of combinatorial simulated annealing technique 
(Press et al. 1992). The obtained configuration of the den- 
sity and velocity distributions should have the least dissi- 
pation and therefore will exist longer than all others. 
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